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Abstract 

We demonstrate that the inversion method can be a very useful tool in providing an infrared stabilization of 
3D gauge theories, in combination with the mass operator A 2 in the Landau gauge. The numerical results will 
be unambiguous, since the corresponding theory is ultraviolet finite in dimensional regularization, making a 
renormalization scale or scheme obsolete. The proposed framework is argued to be gauge invariant, by showing 
that the nonlocal gauge invariant operator A 2 ^, which reduces to A 2 in the Landau gauge, could be treated in 3D, 
in the sense that it is power counting renormalizable in any gauge. As a corollary of our analysis, we are able to 
identify a whole set of powercounting renormalizable nonlocal operators of dimension two. 



1 Introduction 

3D gauge theories are not only of a pure theoretical importance. For example, they naturally appear as the very 
high temperature limit of their AD counterpart H), while 3D QED can be used as an effective theory describing 
high temperature cuprate superconductors |2] . 

When studying 3D gauge theories, a key observation is that the gauge coupling g 2 becomes a dimensionful quantity. 
A positive consequence is the ensuing superrenormalizability of 3D gauge theories, meaning that the ultraviolet 
sector is relatively well behaved. Unfortunately, this also inflicts a serious problem in the infrared [3 |. A simple 
dimensional counting allows us to understand this problem intuitively. At consecutive orders in perturbation theory 
in the massive coupling g 2 , an increasing number of inverse momenta is required to obtain a specific dimension 
of e.g. a particular Green function under study. Consequently, at increasing order of perturbation theory, the low 
momentum region becomes more and more divergent. 

A natural solution to this problem might be the dynamical generation of a mass m, so that a perturbative expansion 

2 

in the dimensionless parameter £- might emerge, ensuring a safe infrared limit. This has stimulated several studies, 
ID |6] H E] [5] to quote only a few. A common feature of the approaches of for example J4j [5] 17) is that a 
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dynamical gluon mass is derived from a certain gap equation, constructed in a particular approximation scheme, 
bearing nontrivial solutions. 

The aim of this article is to investigate 3D gauge theories in the presence of the gauge invariant nonlocal operator 

A^iVTy^^J^f. (1) 

In the first part of this paper, we shall prove its power counting renormalizability, making the mass term »i 2 A 2 nin a 
gauge invariant candidate for an infrared regularization. As a nice byproduct of this proof, we shall be able to iden- 
tify a whole class of gauge invariant nonlocal operators which also enjoy the property of being UV powercounting 
renormalizable. However, we shall discuss why A 2 ^ plays a preferential role, since we must also take into account 
potential infrared problems. In the second part we shall employ the inversion method to get a meaningful pertur- 
bative expansion for 3D gauge theories when the regulating mass coupled to A 2 nin is brought back to zero. This is 
done in the case of the Landau gauge, as A 2 ^ then reduces to the local operator A 2 , which has already been studied 
before, revealing several interesting properties ifTolfTTI . We end with a discussion of the results. 



2 The UV power counting renormalizability of A^ in in 3D 
2.1 Preliminaries 



We intend to use the following action in 3D Euclidean space time 



S = / d X \ ~A F l + bd ^f + Td ^ C + i m A l ' 



(2) 



which corresponds to a Yang-Mills theory plus Landau gauge fixing, supplemented with a regulating mass term. It 
was not only proven that this action is renormalizable to all orders of perturbation theory, but even that it is finite 
in dimensional regularization, i.e. no renormalization is needed [ lOl fTTI . 

There are 2 remaining questions to be answered. The added mass term « A 2 does not appear to be gauge invariant, 
and the hitherto free mass parameter m should be of a dynamical nature. Since the 3D gauge coupling carries a 
dimension, one expects that the dynamics of the theory will dictate m°< g 2 . 

The operator A 2 has a gauge invariant meaning in the Landau gauge, 3A = 0. Indeed, the gauge invariant operator 
(Q3 can be rewritten as a perturbative series lfl2l 



= ~ 2 l^ 
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(3) 



and clearly A 2 ^ — A 2 when 3A = 0. 



Recently, in 4D, it has been argued that the highly nonlocal operator A 2 ^ might be handled in a perturbative 
fashion [ 13 1. The idea is to use the following termwise gauge invariant representation of A 2 ■ fl4l . 
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= Tr / d\v ( FnvjpFfjv + Hg-^F^ 

r i 



iDyFx, 
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added to the action in the following form 



m Tr / d x { F^j^F^ + ligj^Fx^ 



j^DkFkx, j^D v F VI j 



0(F 4 ), 



-Jy2.D K F^X, j^DvFvfj 



(4) 



-2ig 



0(F 4 ), 



or in a more condensed notation 



S„, = m 



! Tr J d\x(o 2 +g0 3 +g 2 4 + ...) , 



(5) 



(6) 
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where the 0# are gauge invariant functionals of the field strength F and co variant derivative D. As noticed in |fl3l , 
this expansion can be seen as one in operators with k legs where k counts the lowest number of gluon legs present 
in the operator Ok- 

Since the series (|6]l or (f3]l contains infinitely many nonlocal terms, it appears to be beyond our possibilities to show 
that such a highly nonlocal operator might be renormalizable. In general, the utmost we could achieve is to show 
the renormalizability up to a certain (low) order. Therefore, we consider the expansion ©■ Each of the nonlocal 
gauge invariant terms could be studied separately. Such an approach was successfully employed in AD in the case 
of the first term, F-^F lfl5l [161 . At the cost of introducing a set of auxiliary bosonic and fermionic fields, the 
action with the nonlocal operator added to it, was cast into a local form. Using the many Ward identities of the 
resulting action, we were able to prove the renormalizability to all orders of perturbation theory. We started from 
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recast into 
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with B,B a pair of complex bosonic antisymmetric tensor fields in the adjoint representation and G,G a pair of 
anticommuting antisymmetric tensor fields, also in the adjoint representation. We succeeded in constructing a 
gauge invariant classical action S c i containing the mass parameter m which is renormalizable. This was proven to 
all orders of perturbation theory in the class of linear covariant gauges, and explicitly checked up to two loop order 
lTT5l[T6l . In particular, this action reads 
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We draw attention to the fact that an additional quartic tensor coupling \ abcd , as well as two new mass couplings 
A-i and A-3 had to be introduced in order to maintain renormalizability. This fact obscures the identification with the 
original operator FjpF, and hence with A^ in . Nevertheless, at one loop, these extra couplings are not yet relevant 
in the renormalization of A^ in , as found in [13]. The renormalizability was explicitly confirmed at 1-loop in a 
general linear covariant gauge, with a gauge parameter independent anomalous dimension. The retrieved value 
did coincide with the already known result for the anomalous dimension of A 2 in the Landau gauge (TTl [T8l , as 
expected from the gauge invariance of A^ in and the perturbative equivalence with A 2 . In this sense, the result of 
|fl3l is very interesting as it provides evidence that A^ in could be consistently used at least at lowest order. Since it 
is gauge invariant, one can opt to work in the Landau gauge, where it reduces to a single local operator, which can 
enter the OPE for example. Needless to say, this also considerably simplifies practical calculations. However, so 
far, the analysis was restricted to lowest order. Beyond the 1-loop approximation, little is known. Things inevitably 
will become complicated since the new coupling \ abcd will explicitly enter the analysifl and evidently, when A^ in 
would be renormalizable, its anomalous dimension is not supposed to contain any new couplings. 

In principle, a completely similar pathway could be followed in the 3D case. One could investigate whether the 
consecutive terms in the expansion (f5]l are renormalizable, by introducing extra fields etc. Almost needless to say, 
this is still a very cumbersome job. For the first term O2, this would amount to a 3D analysis of its localized version 
similar to 4D |fT31[T6l . The Lorentz structure might be simplified a bit by using the dual vector field = jE^ka/ScA, 
and analogs for the localizing fields. A potential caveat would be the emergence of the extra couplings, which 
again make it obscure (or even make it impossible) to say that O2 itself is a renormalizable operator. However, it 
might be very well possible that the massive 3D version of (O without extra couplings (A,j = A-3 = \ abcd = 0) is 



'The 2-loop anomalous dimension of m is X abcl1 dependent 1161 . 
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renormalizable. We recall that these couplings were originally introduced to absorb generated new counterterms. 
At 1-loop, the 3D theory ought to be finite, as the "master integral" is finite in dimensional regularization IfTTl . 
beyond 2-loops no new counterterms can arise (similar arguments as in IfTTl ). so the only possible source of these 
extra couplings would lie at 2-loop order. In principle, this could be checked similarly as done for A 2 in the Landau 
gauge [ 1 1 j. In this work, we shall follow a slightly different route, as the situation might be more appealing in 3D 
due to the superrenormalizability. 

2.2 Inductive proof of the UV power counting renormalizability of A^ in 

In this section, we shall establish the UV power counting renormalizability of A 2 ^, As this proof will turn out to 
be rather technical, let us give a brief sketch of the main argument. We notice that in the expansion (0J, vertices 
will appear with an arbitrary power of the coupling g. Hence, as g becomes dimensionful in 3D, vertices with a 
certain power of g will necessarily induce a certain power of compensating momenta in the denominator of the 
analytical expression corresponding to a Feynman diagram containing this vertex in order to obtain the correct 
dimensionality. Therefore, one may suspect that this A^ in could be UV power counting renormalizable in 3D. 

Let us now put the previous line of intuitive reasoning on a more formal footing. We start the discussion from 
the operators Ok defined in expression ||6). These operators give rise to a set of new vertices. We shall give 
a diagrammatical inductive argument that these vertices are sufficiently suppressed in the UV such that no new 
ultraviolet divergences will appear. The already present countertermfl of the starting action should be sufficient to 
render the complete theory finite. 

The mass dimension of the operator Ok is actually given by dim [OjJ = 2 — |, since we have dim[g] = j. We are 
in 3D, and each vertex is already multiplied by m 2 . Consider a vertex V, with i (i > 2) gluons legs present. We 
are interested in the high momentum influence of such a vertex V,, i.e. we are interested in its UV "cost" for the 
renormalizability analysis. Since dim[A A ,] = j and a vertex with i legs is multiplied by g'~ 2 , we can conclude that 
Vi ~ gPT' where we represent in general (a combination of) momenta by the rather symbolic notation Q. 

Consider now a random renormalizable set of diagrams of the original theory. More precisely, we consider an 
arbitrary set of (connected) n-point functions. Since the original theory is supposed to be renormalizable, its n- 
point functions are finite after including all the necessary counterterms. To begin, we want to attach a single vertex 
Vi onto these n-point functions, in order to obtain a certain n'-point function. One can check that we obtain all 
possibilities just by connecting a series of external gluon legs. As a consequence there are only two possible 
operations we can undertake: 

1 . We can attach an external gluon leg of the vertex V, to an external gluon leg of an arbitrary n-point function. 
Take £ the number of gluon legs of V, attached to original diagrams, then 1 < £ < i. 

2. We can also glue two gluon legs of the vertex V, together. We assume that in total 2s legs are pairwise closed 
on each other to form loops. 

After carrying out these operations, e gluon legs of the vertex Vj remain and will serve as external legs. Hence, we 
have 

2s + £ + e = i. (12) 

We notice that the above procedure will generate all possible n'-point functions of the new theory in which a single 
new vertex Vi has been used, as the original starting set was arbitrary. 

We shall now search for an estimate of the UV "weight" W caused by the new vertex insertion. We shall work in 
a "worst case scenario", i.e. we always consider the case that the UV behaviour is the worst of possible occurring 
scenarios. In order to make things as comprehensible as possible, we shall explore all possible occurring scenarios 
one by one. 

2 As far as these are existing of course, in case the starting theory would be finite. 
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2.2.1 The vertex 



As already explained before, the vertex V, will induce a weight 



1 



Wv ~ — f . (13) 



2.2.2 The s loops 

As it is easily checked, the 2s legs closing on each other will generate the following contribution after integration 

* - d4) 



(G 2 ) s 



2.2.3 The attachment of the £ legs 



Each of the £ legs of Vj shall be glued to an external leg of an original diagram, which is part of a renormalizable 
set. Such an external leg is coming from a 3- or 4-point verted of the original Yang-Mills theory. We will consider 
case by case, whereby every case is determined by the number of legs of V, that "arrives" at the same vertex of the 
original diagram. 

• Case 1: the single attachment 

Let us assume that there are £\ spots at which a single leg arrives. There are 4 scenarios, as shown in Figure 
1: 

1 . £\ a by using a 3-vertex with 2 external legs 

As clearly depicted in Figure [T(a)} this corresponds to using £' la times a tree level 3-point vertex as 
starting point. Each time, one leg of it is glued to a leg of the new vertex V,. The UV weight is obtained 

as 



~ I 7^2 ) ~ — ■ (15) 



The corresponds to the extra propagator caused by gluing 2 legs together, the Q comes from the 
3-point vertexQ, and there are no loop integrations possible in this case. 
£' lb by using a 3-vertex with lorO external legs 

In comparison with the previous case, one of the legs no longer serves as an external leg, but becomes 
connected itself to another vertex of the old diagram. This possibility is depicted in Figure \Ub)\ where 
the grey area stands for any other diagram^. In the current case, one can create additional loops which 
will influence the to-be-derived weight factor. In addition to the foregoing weight factor, we shall also 
encounter extra loop integrals. In particular, 

pi 

w i[t ~ (| d^e) 16 ~ (e 2 )V (i6) 

However, this counting procedure is not as fine as desired to obtain a reasonable final result, being a 
sufficiently suppressed UV weight. Fortunately, there is an alternative way of obtaining such kind of 
vertex insertion, which allows for a refined weight factor. Namely, taking a look at Figure [1(5)1 we see 
that the original diagram can be thought of as another diagram of the old theory that we have cut open, 
added an extra old 3-vertex to it, and then we have attached to this extra vertex the legs coming from 
the new vertex V,-. The upshot of this viewpoint is that now, we can take into account that the "cut & 



3 The 3 refers to the 3-gluon and ghost-gluon vertex. For simplicity, we have omitted the quarks for the moment. 
4 Notice that this is a worst case scenario, as the ghost-gluon vertex does not carry a momentum factor. 

5 In order not to overload the picture, we did not draw (the) other leg(s) connected to the grey blob assuring that the final diagram would be 
IPI. 
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(a) (b) (c) (d) 



Figure 1: The possible Feynman diagram configurations for Case 1. Fat lines refer to the new vertex V;, while normal 
lines to the original theory. 

paste" operation on the new original diagram brings an extra propagator, viz. ^jy, into the game. As 



such, we obtain 



rather than (fT6l l. 



'■lb 



Mi' ~ VQ—^Q ~ (1)^ = 1 (17) 



A word of caution on the emerging loop integrals: it should be understood that these loop integra- 
tions are performed at the end, when all attaching operations are done. However, for the purpose 
of counting the UV weight, we have distributed them over the several subcases, since this does not 
influence the formal power counting and allows for a more efficient bookkeeping. 

3. t'[ a by using a 4-vertex with 3 external legs 

This case (see Figure |l(c)| i can be treated in an analogous fashion as the £' la case, albeit no Q will 
appear since a 4-vertex does not contain momentum dependent factors. We find 

<W f n ~ — —sr. (18) 

4. I'L by using a 4-vertex with 2,1 orO external legs 

This is analogous as the l\ b case, but without the Q, thus yielding 

/ f 1 \ hb 1 

- [I^wy) ~w (19) 

The corresponding diagram is shown in Figure [T(d)| 
Case 2: the double attachment 

Let us assume that there are li spots at which 2 legs arrive (see Figure 2), whereby evidently each time a 
loop is created. Also here, 4 configurations arise: 

1 ■ ^2a by using a 3-vertex with 1 external leg 
The following UV weight is found 

^ ~ {h 3Q w? Q T~ l - (20) 

For the benefit of the reader, let us again explain the origin of the different components in the pre- 
vious weight factor. There is a loop integral, the two propagators building the loop and an extra Q 
corresponding to the 3-vertex (cfr Figure |2(a)] i. 
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(a) (b) (c) (d) 



Figure 2: The possible Feynman diagram configurations for Case 2. 

2. £' 2b by using a 3-vertex with external legs 

Using a slightly adapted argument, namely a pure "paste" one, we obtain in this case (see Figure [2(b)] ) 

Depicting the procedure for the reader: we attached 2 legs of the new vertex V,- to a single external leg 
of the grey blob by means of 3-gluon vertex. 

3. £'2 a by using a 4-vertex with 2 external legs 
Now one recovers (see Figure [2(cJ| > 

~ {h' Q w?T^w (22) 

4. ^2fc fl 4-vertex with 1 or external legs 

For the fourth double attachment scenario (see Figure [2(d)] ), we have 

X f" 



which is obtained by the "cut & paste" logic. 



• Case 3: the triple attachment 

Let us assume that there are It, spots at which 3 legs arrive. In this case, we observe 3 options: 

1 . £' 3 by using a 3-vertex 

This is the simplest case, as no external legs are available. We find 

1 \4 



*4 ~ [jmy-^Qj ~q*i. (24) 

There is a double loop integral, as it immediately follows from the diagram displayed in Figure [3(a)| 

2. £" a by using a 4-vertex with 1 external leg 
From Figure [3(b)] we find 



1 



(Q 2 ) 



r" 



3. £'^ b by using a 4-vertex with external legs 



Once more employing the "cut & paste" argument leads us to 

1 



3/) 



The corresponding diagram is shown in Figure [3(cjl 
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Q 




(a) (b) (c) 

Figure 3: The possible Feynman diagram configurations for Case 3. 



Figure 4: The possible Feynman diagram configuration for Case 4. 

• Case 4: the quadruple attachment 

Finally, we can assume that there are £4 spots at which 4 legs arrive. Evidently, only one possibility pops up: 

- (/( d3 £) 3 (^) 4 ~e 4 . (27) 

This was obtained analogously as in Case 2.2. 
For later use, we mention that 

I = 4 + 4* + 4, + 4* +2 (4, + 4* +4, + ^) +3 (4 + 4^ (28) 

=h =h =h 

which is readily checked. 

We are now ready to combine all the obtained information into a single estimate for the UV weight, 

W — ^i<Wv<W s <W f i Wei Wfi 'Win W f i <W f i Win Wen Wfi Wen Wen W(, . (29) 

(J-> lfl lb la lb '2a 2b '2a '2b 3a '3a 3b 4 

We introduced a factor ^j, which serves as a "correcting" weight factor related to local momentum conservation 
at the vertex V,. During the derivation of the different weight factors, we have always assumed that the introduced 
loop integrals were independent, but momentum conservation at V, will at least kill one of these. Let us also 
mention here that the e external momenta of V, will get related to the other external momenta of the final diagram 
by means of global momentum conservation. 



8 



Simplifying expression ( |29l , we find for the total weight: 



with k = (i - 2) + 3 -s + 4, + 2^ + 4, - 4 6 + 4, - 4 - 4ft - ^4 • 
The power K can be further simplified by means of (fT2b and d28l ). which leads to 

k = 1 + s + e + 2/ la + e xb + W[ a + 21% + 24 + 4, + 311 + 24, + 2/3' + 3£ 3 ' a + 2^ + 31, > . (31) 

Recapitulating, the power of g is strictly positive, meaning that the UV weight <W , (f30b , of the new vertex is under 
control very well when Q ~ 00. 

So far, we have ignored the presence of quarks. The counting analysis will however remain valid even in that case. 
The only cases in which the 3-point quark-gluon vertex will appear and influence the power counting, correspond 
to the diagrams displayed in the figures [T(a)| or |l(b)| However, the situation corresponding to the case |l(a)| is even 
further improved since the quark-gluon vertex does not contain an explicit momentum dependence. For the case 
of I l(b)| we notice that the extra quark propagator will behave as ^ instead of ^j, but since there is no Q coming 

from the vertex itself, the eventual counting remains unaltered since "^z<2 = g"- 

To close the inductive argument, we can of course repeat the previous argument when we allow a second, third, . . . 
new vertex of the type V, into where we are now: the theory plus one single vertex. The latter one has just been 
shown to behave well in the UV. We conclude that any number of new vertices added to the theory will only induce 
UV harmless additional (pieces of) diagrams. 

Before closing this section, a few more words are to be devoted to the i = 2 case. Setting i — 2 (and thus s = e = 0) 
in the UV estimate d30l ), we see that also the 2-legged insertion is UV safe when used to construct new diagrams 
starting from an originally renormalizable theory. In this case, this would be massless 3D YM, whose perturbation 
theory is unfortunately ill-defined. We can however circumvent this problem. In every "conventional" gauge, there 
exists a renormalizable (gauge variant) mass operator: the Landau [ 19 1, the linear covariant EDI , the Curci-Ferrari 
ED . the maximal Abelian [21 1 and a class of nonlinear covariant gauges [22|. Strictly speaking, this was proven 
only in AD, with the exception of the Landau and Curci-Ferrari gauges which were also explicitly analyzed in 3D 
11101 . but the algebraic renormalizability analysis of the mass operators does not really depend on the space time 
dimension, and all the relevant Ward identities will remain valid in 3D. Continuing our reasoning, we choose a 
gauge to work in and add an IR regulating mass term to it. Starting from this theory, we can apply the foregoing 
arguments of this subsection also for the 2-legged insertion, and consequently conclude that A^ in is power counting 
renormalizable in any of these gauges. Once this is established, we can drop again the temporarily introduced gauge 
variant mass term and immediately add the gauge invariant mass term °c A^ in to the action. 

Summarizing, we have thus demonstrated that A^ in should be power counting renormalizable in 3D. The gauge 
invariant mass operator A^ in should thus be consistent with the ultraviolet renormalizability. For practical calcula- 
tions, it would still be technically challenging to calculate with A^ in even in the Feynman gauge. However, since 
it is explicitly gauge invariant, we do not make any sacrifices choosing to work in the clearly preferable Landau 
gauge, in which case we have the relatively simple action (0 at our disposal. In this case, we also do not have 
to worry about the overabundance of ^7 in the expressions © or ©, which could cause IR troubles during the 
calculation of the Feynman integrals. Due to the gauge invariance and the subsequent choice of the Landau gauge 
where no potentially dangerous 4i terms are present, it is obvious that these possible infrared divergences cannot 
influence physical quantities (they should cancel out if occurring in intermediate results). We shall come back to 
this issue at the end of section 3. 



2.3 A more general class of powercounting renormalizable gauge invariant 2D operators 
and the preferred role of A^ n 

As the attentive reader might have noticed, we can extract more interesting information from our proof than simply 
the UV powercounting renomalizability of A^ in . In fact, all the operators O, are separately UV powercounting 
renomalizable. A fortiori, so is any (infinite) linear combination of those. However, there is more to the story in 
3D than just UV powercounting renormalizability, we should also take the infrared safety into account. First of all, 
O2 should be present as this is the only one that will give rise to a mass in the gluon propagator, which can serve 
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as a natural infrared cut-off. Secondly, we should also be aware of the potential infrared danger caused by the ^-'s 
in the interaction terms. It is exactly our point that by using a particular series of these operator, viz. A^ in , one can 
motivate that no infrared dangerous terms will occur when calculating gauge invariant quantities, as we have done 
at the end of the previous subsection. 



3 Removal of the regulating mass parameter m 
3.1 The inversion method 

Albeit that we have regularized the 3D gauge theory in a gauge invariant fashion by the introduction of the mass 
m, it is still a mass introduced by hand. The next goal is to get rid of this arbitrary parameter. 

We shall explain the inversion method with the example of the pole mass. Consider the one loop gauge boson 
polarization tensor 7t^, which can be decomposed in the traditional way 

<(, 2 ) = 8«* ( 8, v - m) n(q 2 ) + 5^co(, 2 ) . (32 ) 



) — « I «//v — ) n Vi > u — T 

\ *i J y[ 

It is then easily shown that the corresponding shift in the tree level mass will be given by 

m 2 -^m 2 + %(q 2 ). (33) 

The pole mass (squared) m 2 is defined as that value of (— q 2 ) such that q 2 + m 2 + n(q 2 ) = 0. As we are working 
with a perturbative series, this can be solved in an iterative way, so that at lowest order 

m 2 p = m 2 + g 2 %i (-to 2 ) , (34) 
where Tti (q 2 ) is the one loop contribution to the self energy. 

Having found an expression for the pole mass m p in terms of the regulating mass m, 

r 

m m 



m 2 =m 2 ( 1+fli— +a 2 ^r + ... ) =m 2 A{m 2 ), (35) 



the question remains what we must do with it, since m should in principle become zero again at the end of any 
calculation, to restore equivalence with the original starting massless YM action. One option is to simply set m = 
in ( |35T >. However, as outlined in 11231 , (certain) nonperturbative effects can be taken into account when the series 
( |35T > is inverted as follows, 

m 2 = m 2 ( 1 + bi — + Z? 2 ^2 + • • • J = ( w p) • ( 36 > 

Apparently, m — can now be realized in 2 ways: by setting m p — corresponding to the perturbative (potentially 
ill-defined) solution, but also by solving the gap equation 

SK)=0, (37) 

which can give rise to a nontrivial solution m 2 ^ 0, whereby that nevertheless m = 0! 

We used the example of the pole mass m p to explain the inversion philosophy, but also other quantities Q, could 
be handled: the regulating mass m is introduced to ensure a meaningful perturbative series for Q,(m), and after 
inversion a meaningful (finite) result can be found for Q, even for m = 0, by solving the gap equation m(d) =0. 
Hence, it appears that the inversion method can be a very useful tool to obtain results in superrenormalizable 
quantum field theories, which are plagued by infrared instabilities. 

3.2 Explicit calculations 

We now turn to an explicit computation. We have calculated the gauge boson self energy at one loop by evaluating 
the contributing four Feynman diagrams in three dimensions where we also include massless quarks. The diagrams 
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are generated by the Qgraf package, [24], and converted to FORM input notation where Form is a symbolic 
manipulation language, ||25ll . The self energy is then reduced to a set of master one loop integrals which have 
previously been determined in [26|. Whilst these have been deduced in dimensional regularization in d = 3 — 2e, 
the relevant one loop integrals are in fact finite in three dimensions. Hence we find that the transverse part of the 
self energy is 



2 ) = 



1 Jq 2 

q 2 Arc Tan I - — — 

2 m 




i^ArcTanf^U-L 
8 m z \ m J 32 



for the transverse component, and 

U(q 2 ) 



1 

128 



-,2\2 









{-re- 


1 

~8 V 


1 

1 2 
q L 




1 

76 



(q 2 ) 2 



ArcTan 



(38) 



(39) 



for the longitudinal component. The reader might be a little confused, as the computed 1-loop self energy Ttf^ is 



apparently not transverse. However, the Ward identity in the massive Landau case does not predict a transverse 
self energ}@. Details can be found in the Appendix. 



3.3 Back to the pole mass 

Let us return to the determination of the pole mass, at 1-loop given by 

m 2 = m 2 + g 2 %i(-m 2 ) . (40) 

A few complications arise. A first one is the appearance o£| ArcTan(/) = ;°° in 7c(— m 2 ), but fortunately, these 
terms cancel amongst each othei0. Less fortunately, the presence of \/? results in a complex valued pole mass. 
We find 

.2 _. _2 *'^.„2.„ N f 3 63 \ , ;Ar 1 , 



p = ~ l T F N f8 m +~ ^ - I^ ln3 J r« + W— g^mH- . . . , (41) 

since ArcTan (=) = ^ln3. Performing the inversion and solving the induced gap equation, we obtain 

i 9 N f 3 63 \ o 1 o 
»V - - g T F N fg > + -{ T2 -^) 8 2 + iN—/ + .... (42) 

It is interesting that, using a different approach and different mass operators, the works ||6]|7) also report a complex 
pole mass in some cases, depending on the employed (gauge invariant) mass operator. 



Evidently, the propagator (connected 2-point function) is still transverse. 



6 

7 We recall that ArcTan(z) = i In |±|. 
8 As a matter of fact, also the longitudinal part is finite at q 2 = —m 2 . 
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We recall here that 3D gauge theories are also confining. This means that the gluon itself is not a physical particle, 
and as such, unitarity should not be expected at the level of the elementary gluon excitations. Therefore, we 
certainly do not claim that we have obtained massive gauge bosons with 3 physical polarizations, even if the 
eventual pole mass would have been real-valued. Notice that unitarity at the level of the gluons was also not 
mentioned in e.g. (6] [7J. The main point of these papers and the current work is to find a way to regulate the 
3D gauge theory to allow for a consistent expansion. Once this is done, one could try to have a look at the 
physical excitations, which are supposed to be massive glueball states, a fact supported by the 3D lattice data 
11271 . For example, one could try to construct a bound state of gluons using standard techniques, or one could 
study gauge invariant correlators like (F^ v (x)F^ v (y)), since the operator F 2 V has the correct quantum numbers to 
create/annihilate a scalar glueball. This is however beyond the scope of this paper, but let us only underline that an 
explicit calculation of (F 2 V {x)F* (y)) would also be plagued by infrared singularities in 3D, unless some regulating 
mechanism is being provided. 



3.4 Applying the inversion method to the gluon propagator 

As a second example, we shall now determine an estimate for the gluon propagator. In our conventions, the gluon 
propagator reads 

D ^ = 2^ 2lW 2 2^ ' < 43 > 

which we have calculated explicitly to 1-loop order. For any value of the Euclidean momentum p = p*, we can 
invert (l43l to get m 2 = f (Z)(p*)), and then solve J (D(p*)) = in order to obtain a value for D(p»). Since we 
have 

= p 2 + m 2 + g 2 %i(p 2 ,m 2 ) (44) 



D{p) 

at lowest order, where the value of Jti (p 2 ,m 2 ) can be extracted from ( UoT l. we can solve for m 2 in an iterative way 
as follows 

m 2 =^p)-g 2 Ti l (p 2 £), (45) 
where we put !;(/>) = j^-t — p 2 for notational simplicity. Assuming that ^* (p) obeys 

Up)-8 2 Mp 2 Mp))=0, (46) 
we arrive at the following 1-loop estimate for the 3D gluon propagator 

D{p) = 2 * - f ; ■ (47) 
P z + %*ip) 

We plotted Dip) as a function of p in Figure 5. 

Since we are still using a series expansion, one might wonder if we have any control over this expansion? We recall 
that the formal counting of the orders in the expansion is done by g 2 , which unfortunately carries a dimension, 
hence g 2 is not really suitable as expansion parameter. However, taking a look at the inverted series d46l ). we notice 

that we can say that -A= will emerge as a natural dimensionless expansion parameter at each ordeiQ. In addition, 

\/%* 

in a D-dimensional space time and in the absence of quarks (Nf = 0), the coupling constant is always accompanied 

N 
(4tc)°/ 2 



by an additional suppressing factor jr^m due to the loop integrations. In Figure 6, we have therefore plotted the 

quantity y(p) = f —. in function of the momentum p. We recognize this is a rude way of estimating the 

acceptability of a perturbative approach, but at least we can be satisfied that y is sufficiently small if we do not 
come too close to zero momentum. 



3.5 A few extra words on A^ in beyond the Landau gauge 

We recall that A 2 is not gauge invariant, but in the Landau gauge it equals A^ in , and we motivated that the latter 
gauge invariant quantity is power counting renormalizable in 3D, assuring that we in principle considered a gauge 

9 This amounts to a dynamically realized version of the naively expected ^ in the presence of a regulating mass parameter m. 
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invariant regularization. We draw attention to the fact that the ghost in the Landau gauge is still massless, and might 
introduce additional infrared instabilities^- In order to overcome these, one might consider using the Curci-Ferrari 
gauge, a generalization of the Landau gauge, in which case the ghosts also attain a mass ifTOlfTTl . At the end of 
such a calculation, the limit to the Landau gauge can be considered. However, since the ghost mass is explicitly 
gauge parameter dependent, it will not enter gauge invariant quantities, and as such, the Landau gauge could be 
immediately used to calculate gauge invariant quantities without risking extra infrared divergences. 

The question remains what to do if we would like to calculate e.g. the gluon propagator D(p) in the linear covariant 
gauges with A^ in ? Infrared singularities coming from the 4r's in (0 can be avoided by using 



S m = m Tr / d x I F, 



//V 



1 1 



-2ig 



D 2 + p 2 



F-Kfj 



1 1 

jD k Fks,, , _2^MV 



D 2 + p 2 K KV 'Z) 2 + p 2 



D2 + p2~™D2 + p2~ V ' V " 

0(F 4 ) 



(48) 



instead of ©. We introduced an extra mass parameter p 2 , but neither the gauge invariance nor good UV behaviour 
are compromised by this[[]} The potentially dangerous ^ will get replaced by the IR safe ^rj^i ■ As such, we will 

get a well defined D(p,m 2 ,p 2 ). Once this is done, one can check whether the limit p 2 — * exists or not. In the 
Landau gauge, this should be the case as already explained before, in which case we are back to A 2 -;- . If not, one 
can perform a first inversion with respect to p 2 to find a sensible D(p,m 2 ), starting from which a second inversion 
can be done to ensure m 2 = 0. 



10 These were not observed during the two-loop calculations reported in ill), despite the presence of massless ghosts in the Landau gauge. 
1 1 Notice that a similar kind of modification would not be possible using the series {3}, as this would spoil the gauge invariance. 
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4 Discussion 



The purpose of this paper was to illustrate that inclusion of the operator A 2 in combination with the inversion 
method allows for a consistent (i.e. infrared protected) perturbative expansion of basically any quantity. However, 
there are other major sources of effects beyond (regularized) perturbation theory. For example, we can mention 
the existence of Gribov (gauge) copies in the Landau gauge. Trying to restrict the integration measure in order 
to take these copies into account, can have a profound influence on the gluon propagator, analogously as in AD 
ll28l . The restriction introduces another mass scale y in the theory [29|, and we can expect that y °< g 2 in 3D. 

A standard Gribov-like propagator looks like p <t + ^, , clearly having 2 complex poles at p 2 = ±/y 2 . The Gribov 
restriction is important as a theoretical tool to find a violation of positivity in the gluon propagator, indicative of 
confinement ll29l . In this light, dynamical gluon mass scales, complex or real, should not be directly related to 
massive "physical" gluons, as these are confined and hence unphysical. Therefore, a complex pole mass as found 
here or in other works EG) is not necessarily a catastrophe. We recall that 3D gauge theories, although that the 
classical/perturbative interaction potential is already mildly (viz. logarithmically) confining, also display "true" 
confinement trough a linear potential, see e.g. [27|. The origin of this piece of the potential is not evident. 

When we compare our propagator displayed in Figure 5 with the lattice result of ll30l . then it is immediately clear 
that the big difference is located in the deep infrared: the lattice results indicate a. finite gluon propagator near zero 
momentum. But at larger momenta, the inversion mechanism described in this paper which regulates the theory 
can give acceptable results. 

To conclude, it would be recommendable to pursue e.g. an analytical study in 3D based on OTI taking into account 
the existence of Gribov copies, and find out whether a more quantitative agreement with the available lattice data 
can be found also in 3D \29, 32]. At the same time, it can be investigated whether the Gribov scale y would serve 
as a natural IR regulator. 
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A The Ward identity for the gluon self energy 

We start from the complete classical action 

E = S YM + S Landau + Sext + s j d 3 * 



J d 3 x (\f^ v + b a \Al + <fd M D?c b + l -JA"^ + TA; V - Slfflf c b + 8 -f bc L u c b A , (49) 



supplemented with the necessary extra (external) source terms, e.g. J which is used to couple the operator A 2 to 
the theory. Moreover, s denotes the usual BRST symmetry generator, 

sAl = -(Dp) a , sc a = l -gf abc (*c c , sc" = b\ sb a = 0, (50) 
extended to the sources by means of 

sx=J, sJ=0, sQ. a M = Q, sL"=Q. (51) 



The corresponding Slavnov-Taylor identity reads 
5(E) = Jd 3 x 



3 „/5E 5E , 5E5E +/7 „5E +y p 



\8Q»8A<* hL"hc a 6c" 8t 
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The Ward identity for the vacuum polarization can be derived from this Slavnov-Taylor identity, suitably extended 
to the quantum level. At the one loop level, one has 



r = E+fir' 

so that the 1-loop Slavnov-Taylor identity becomes 



(53) 



/A 



3 /sr'SE se sr 1 sr 1 se se sr 1 sr 1 sr l \ 

3 ~ ' 1 1 1 1 b"—+J— =o. 



5£2£ 8A« 8£2£ 8A£ hL a 8c" 8L" 8c" 8c" 8x 



Usini 



,12 



8T^ 



i 8T^ 

' 8Z« 



(1/ 



«&c c 6 c c ) . r 



i 8T^ 



[(W) -r]' 



(54) 



(55) 



we derive 



- (Dfc 6 ) 1 



i 8E 



8A? V " SAf, LV2 7 



i 8E 

8c" 



f/" / ^0£+^(W')-r] 1 + ^£) °- 



8c" 



(56) 



s2 

Applying the test operator " ^ and setting all fields and sources equal to zero at the end of the operation, 
except for J — m 2 , one obtains the Ward identity for the vacuum polarization nfy 



8 2 r l 



8A«(x)SA b v (y) 



8 2 



Sc a (x)?>A°(y) 



V 



/ 



1 1 



\ 



8c°(x) 



\ 



(57) 



/ 



The first term is trivial, by using partial integration and the fact that we are considering the Landau gauge dA = 0. 
Hence, the vacuum polarization in the massive Landau case is not transverse, but rather subject to the following 
Ward identity at one loop 

/ r / \ v\ 



8c°(x) 



(58) 
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